Jeans analysis of self-gravitating systems in f(R) gravity 
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Dynamics and collapse of collisionless self-gravitating systems is described by the coupled collision- 
less Boltzmann and Poisson equations derived from f(R) gravity in the weak field approximation. 
Specifically, we describe a system at equilibrium by a time-independent distribution function fo(x, v) 
and two potentials <&o(x) and 9o(x) solutions of the modified Poisson and collisionless Boltzmann 
equations. Considering a small perturbation from the equilibrium and linearizing the field equations, 
it can be obtained a dispersion relation. A dispersion equation is achieved for neutral dust-particle 
systems where a generalized Jeans wave-number is obtained. This analysis gives rise to unstable 
modes not present in the standard Jeans analysis (derived assuming Newtonian gravity as weak filed 
limit of f(R) = R). In this perspective, we discuss several self-gravitating astrophysical systems 
whose dynamics could be fully addressed in the framework of /(i?)-gravity. 

PACS numbers: 04.50.Kd, 04.25.Nx, 04.40.Nr 



> 

r-" 
o 



X 



I. INTRODUCTION 



One of the fundamental goals of modern cosmology is 
to probe Einstein's General Relativity (GR) at any scale 
beyond the classical tests that confirmed such a theory 
in the weak field limit and at Solar System level. GR 
is assumed as the standard theory of gravity describing 
astrophysical structures up to the whole observed Uni- 
verse; however there are some inconsistencies at ultravi- 
olet scales (e.g. the initial singularity, the quantum grav- 
ity issue) and infrared scales (e.g. cosmic acceleration, 
concordance problem, flatness problem, galaxy rotation 
curves, large scale structure, massive stars formation) 
that strongly suggest that Einstein's approach should be 
revised or at least extended. Furthermore, astrophysical 
observations of the last decades suggest that new (dark) 
ingredients are necessary to achieve a self-consistent cos- 
mological model. In particular, the observations suggest 
the Hubble flow is currently accelerating, and the sim- 
plest way to explain the cosmic acceleration is to insert a 
cosmological constant (A) in the Friedmann- Robertson- 
Walker Cosmology [lH3|, representing about 70% of the 
total amount of energy. On the other hand, the galaxy 
rotation curves and the large scale structure could be dy- 
namically addressed by introducing huge amounts of dark 
matter (about the 25% of the total matter). Only 5% of 
the cosmic budget is constituted by standard matter as 
stars, neutrinos, radiation, heavy elements and free cos- 
mological hydrogen ad helium. Alternative approaches 
to GR could be pursued with the aim to explain the ob- 
served acceleration and missing matter without introduc- 
ing new ingredients up to now not observed at fundamen- 
tal scales. The so-called /(i?)-gravity is considered as a 
possible, straightforward mechanism to explain the cos- 



mic acceleration without inserting unknown elements as 
dark energy and dark matter, but extending the geomet- 
ric part of the field equations by relaxing the strict hy- 
pothesis that the gravitational action has to be restricted 
to f(R) = R as in the Hilbert-Einstein one (3-l8l. |30|. 

These theories have been investigated both at cosmo- 
logical scales and in the weak field limit [IoIh141 ]. It has 
been shown that a late accelerating behavior can be easily 
recovered and it can be related to an early inflation- 
ary expansion [l6| . Furthermore, modifying the gravity 
action by assuming non-linear Lagrangians, one obtain 
corrections to the gravitational potential which can be 
useful for astrophysical phenomenology at galactic scales. 
In particular, without the introduction of dark matter, 
the rotation curves of spiral galaxies and the haloes of 
galactic clusters can be dynamically addressed (T7rj2ll|. 
Several of these extended models reproduce Solar System 
tests so they are not in conflict with GR experimental re- 
sults but simply extend them [22T - |24l |. 

It is important to stress that f(R) gravity has inter- 
esting applications also in stellar astrophysics and could 
contribute to solve several puzzles related to observed 
peculiar objects (e.g. magnetars, stars in the instabil- 
ity strips, protostars, etc. [34|, HH), structure and star 
formation [32l, l33l|. 

Here we analyze the Jeans instability for self- 
gravitating systems in /(ii)-gravity coupled with perfect- 
fluid matter. The aim is to show that several self- 
gravitating systems, in particular those involved in star 
formation (e.g. large molecular clouds or Bok globules), 
can be exactly addressed in this framework by consid- 
ering the corrections to the Newtonian potential coming 
out from /(i?)-gravity. This fact could constitute a re- 
markable signature to retain or rule out these theories at 
astrophysical level. 
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The paper is organized as follows. In Section [XT] 
the classical theory of gravitational collapse for dust- 
dominated systems is summarized. In Section IIII| we 
discuss the weak field limit of /(i?)-gravity obtaining cor- 
rections to the standard Newtonian potential that can 
be figured out as two Newtonian potentials contributing 
to the dynamics. In Section IIVI we recover the disper- 
sion relation and Jeans mass limit while, in Section |V1 
some self-gravitating dust system are considered in this 
approach. The difference between GR and /(i?)-gravity 
are put in evidence, in particular the Jeans mass profiles 
with respect to the temperature. We report a catalogue 
of observed molecular clouds in order to compare the clas- 
sical Jeans mass with the f(R) one. Finally, in Section 
VH results are discussed. 



DUST-DOMINATED SELF-GRAVITATING 
SYSTEMS 



II. 



The collapse of self-gravitational collisionless systems 
can be dealt with the introduction of coupled collisionless 
Boltzmann and Poisson equations (for details, see J37|): 



df(r,v,t) 
dt 



+ (v-Wr )f(f,V,t) 



V 2 $(r,t) =4ttG / f(r,v,t)dv, 



(1) 



(2) 



where v and r mean 3-dimensional vectors in the spatial 
manifold. 

A self-gravitating system at equilibrium is described 
by a time-independent distribution function fo(x, v) and 
a potential &o(x) that are solutions of Eq.([T|) and (JSJ. 
Considering a small perturbation to this equilibrium: 



/(r, V, t) = f (r, v) + e/i (r, v, t) , 
$(r,i) = $ (r) + e$ 1 (r,t), 



(3) 
(4) 



where e <C 1 and by substituting in Eq. ([T]) and ^ 
and by linearizing, one obtains: 

dfi(f,v,t) _ dfi(r,v,t) 
+ v ■ 



dt 



Of 



. dfnf'r.v) => , , dfi(f,v,t) 



dv 



(5) 
(6) 



fo(v), and the so-called Jeans "swindle" to set $o = 0. 
In Fourier components, Eqs.© and © become: 



-iufi + v ■ [ikfi 
-fc 2 $i = 47rG / hdv 



•7F= ' 
ov 



(7) 

(8) 



By combining these equations, the dispersion relation 

Z dfo 



1 



4nG 



dv 



v ■ k — uj 



■dv = 0: 



(9) 



is obtained. In the case of stellar systems, by assuming 
a Maxwellian distribution function for /o, we have 



fo = 



(2ira 2 



v 

2^2 



(10) 



imposing that k = (k, 0,0) and substituting in Eq.((9]), 
one gets: 



1 _2y/2jjGpo C v x e 2a 2 ^ 



0. 



fccr 3 J kv x — uj 
By setting uj = 0, the limit for instability is obtained: 

,2/ n x ^Gp 2 
k {uj = 0) = 2— = k T , 



(11) 



:i2) 



by which it is possible to define the Jeans mass (Mj) as 
the mass originally contained within a sphere of diameter 
Aj: 



4tt / 1 



where 



Gpo 



(13) 



(14) 



is the Jeans length. Substituting Eq. (fl"4")l into Eq. (|T3^) . 
we recover 



Mj 




(15) 



All perturbations with wavelengths A > Xj are unsta- 
ble in the stellar system. In order to evaluate the integral 
in the dispersion relation, we have to study the singolar- 
ity at uj = kv x . To this end, it is useful to write the 
dispersion relation as 



i - -^wa/3) = o, 



(16) 



defining 



Since the equilibrium state is assumed to be homo- 
geneous and time- independent, one can set fo(x,v,t) — 



W(P) 



'2n 



x-/3 



-dx, 



(17) 
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where ft 



Re 



W 



(-) 

\ka) 



kcr 



and x 



a 



We set also oj = iojj and 



0, because we are interested in the un- 
stable modes. These modes appear when the imaginary 
part of oj is greater than zero and in this case the inte- 
gral in the dispersion relation can be resolved just with 
previous prescriptions. 

In order to study unstable modes (for details, see Ap- 
pendix B in H3|) we replace the following identities 



I 



9 _ 

x e 



rdx 



1 



irfteP [1- erf/3] 



erfz 



dt. 



into the dispersion relation obtaining: 



1 - 



V^fccr' 



X'V2 



1 - erf I 



\y/2kc 



(18) 



This is the standard dispersion relation describing the 
criterion to collapse for infinite homogeneous fluid and 
stellar systems [37| . 



with the trace equation 

3Df(R) + f'(R)R - 2f(R) = XT. 



(21) 



Here, T„ 



-1 8(J=gC m ) 



is the the energy-momentum 
tensor of matter, while T = T a a is the trace, □ = ;cr ;cr 

and f'(R) = The signature is (- + ++) H| ). 

For our purposes, we have to start by setting the right 
approximation in the metric tensor g^ v [3!|: 



-(l + 2$(f,x)) + 0(4) 0(3) 
9n» ~ | | , (22) 



0(3) 



0(2), 



where 0(n) (with n = integer) denotes the order of the 
expansion. It is worth stressing that the expansion pa- 
rameter is c -1 and, in the Newtonian limit, we are as- 
suming perturbations up to c~ 2 . This means that in the 
above expression (|2"2")h we can discard terms of order 0(3) 
and 0(4) that have to be considered in further perturba- 
tion post-Newtonian limit (see Q and references therein) . 
The set of coordinates^ adopted is 2+ = (t, x 1 , x 2 , x 3 ). 
The Ricci scalar becomes 



III. NEWTONIAN LIMIT OF /(i?)-GRAVITY 

As discussed in the Introduction, /(ii)-gravity is a 
straightforward extension of GR by which it is possible, 
in principle, to recover good results of GR without impos- 
ing a priori the form of gravitational Lagrangian, chosen 
to be f(R) — R by Hilbert and Einstein. This means that 
we do not impose a priori the gravitational action but 
it can be, in principle, re-constructed by generic curva- 
ture invariants and then matched with observations (the 
simplest choice in this sense is to take into account an 
analytic function of the Ricci scalar R [i|). However, 
from a genuine mathematical viewpoint, the initial value 
problem of such theories has to be carefully addressed in 
order to achieve self-consistent results (see for example 

M)- 

Let us start with a general class of higher order theories 
given by the action 

A = J d*xy/=g\f(R) + XC m ] , (19) 

where f(R) is an analytic function of curvature invariant 
R and x = is the- usual coupling of gravitational field 
equations Q . The term C m is the minimally coupled or- 
dinary matter contribution. In the metric approach, the 
field equations are obtained by varying (fl9|) with respect 
to g^. We get: 

f'(R)R^ - \f{R)g^ - V„V„/ , (ie)+ 

+g flu Df(R) = XT flu , 

(20) 



R ~ i? (2) (tx) + 0(4). (23) 
The n th derivative of Ricci function can be developed as 

f n (R) ~ /™(i? (2) +0(4)) ~ 

~ / n (o) + .r +1 (o)i? (2) + 0(4), 

(24) 

here R^ denotes a quantity of order 0(n). It is worth 
stressing that the symbol f n (R) means the n th deriva- 
tive of the analytic function f(R). In the following, 
we are going to use the numbers /™(0) of the Taylor 
series. From lowest order of field equations (|20[) . we 
have /(0) = which trivially follows from the above 
assumption (f2"2"j) that the space-time is asymptotically 
Minkowskian. Eas.([2"0]) and (21]) at 0(2)-order (Newto- 
nian level) become: 



R { u ] ~ \- ~ f"(0)V 2 RW = X T t f , (25) 
-3/" (0)V 2 i?^ - R^ = X r<°> , (26) 



1 Here we shall adopt the convention c = 1. The convention for 
Ricci's tensor is R^u = R" p,av while for the Riemann tensor 
is R a pi lv = + .... The affinities are the usual Christoffel 

symbols of the metric = |g'" T (g QCTi/ 3 + g Pa>a - g a p, a ). 

2 The Greek index runs between and 3; the Latin index between 
1 and 3. 
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where V is the Laplacian in the flat space, R). t ' = 
V 2 <I>(i, x) and for the sake of simplicity, we set /'(0) = 1. 
We recall that the energy-momentum tensor for a perfect 
fluid is the 



Tfu, = {e + pju^Uv -pg^v, (27) 

where p is the pressure and e is the energy density. If we 
consider a p erfect fluid of dust (p = 0), we have Rq = 
_ V 2 _goo [30| . Then we have 

r(2) 

V 2 $ - ^- - /"(0)V 2 i? (2) = Xp , (28) 
-3/"(0)V 2 i? (2) - R (2) = Xp, (29) 

where p is the mass densitjB For /"(0) = 0, the stan- 
dard Poisson equation V 2 $ = AnGp is recovered. 

The solution for the gravitational potential $ has 
a Yukawa-like behavior depending on a characteristic 
length. Then as it is evident, the Gauss theorem is not 
valid since the force law is not oc |x| -2 . The equiva- 
lence between a spherically symmetric distribution and 
point-like distribution is not valid and how the matter is 
distributed in the space is very important[___. 

Besides the Birkhoff Theorem at Newtonian level is 
modified: the solution can be only factorized with 
a space-depending function and an arbitrary time- 
depending function. Furthermore, the correction to the 
gravitational potential is depending on the first two 
derivatives of f(R) in R = 0. So different analytical 
models, up to the third derivative, admit the same New- 
tonian general solution. 

Field equations (f2"5| and (f2!J|) give rise to the modified 
Poisson equations for /(i?)-gravity. We know that 

f?(2) ~ 1V72 J 2 ) 1V72 J 2 ) (o n \ 

R y - 2 v 5oo - 2 v 9u ■ (aUJ 

Inserting in the above result the g^ approximations (|22p 
we obtain 

R^ ~ V 2 ($ - *) , (31) 
where ^ is the further gravitational potential related to 

f 21 _____ 

the metric component g u . Substituting in Eqs. 
and (|29|) . we obtain 

V 2 $ + V 2 * - 2/"(0)V 4 $ + 2/"(0)V 4 1' = 2Xp (32) 
V 2 $ - V 2 * + 3/"(0)V 4 $ - 3/"(0)V 4 1' = -Xp (33) 



3 We remember that e = pc 2 . 

4 However, we have to see that being the Yukawa correction a 
decreasing exponential function, the Gauss theorem is asymp- 
totically recovered. In any case, conservation laws are always 
preserved since the Bianchi identities hold. 



By eliminating the higher-order terms, the standard Pois- 
son equation is recovered. Our task is to check how the 
Jeans instability occurs in /(i?)- gravity. 

An important consideration is in order at this point. 
As we pointed out above, we are supposing that the 
space-time is asymptotically Minkowski. However, this 
is against the general idea of /(-R)-gravity which should 
mimic dark energy behavior. This means that the space- 
time should be asymptotically de Sitter. So, in general, 
it is necessary that the f{R) function is expandable at 
R = 0, or even if it is, the interesting asymptotic space is 
nevertheless not Minkowskian (R — 0) but R ^ 0. This 
fact is connected with the assumption that the energy- 
density p is homogenous and asymptotically constant in 
order to leads to de Sitter space-time. This is, at the very 
end, why the so called "Jeans swindle" is needed in New- 
tonian theory. In the present case, p is explicitly written 
in Eqs. (|2^|) and ([29)1 and has to converge asymptotically 
to zero in order to restore the asymptotic Minkowskian 
behavior. In other words, the possible gravitational ac- 
tions have to be chosen so that the condition /(0) = 
holds. 



IV. JEANS CRITERION FOR GRAVITATIONAL 
INSTABILITY IN f(R) GRAVITY 

Our task is now to study the Jeans instability in the 
framework of f(R) gravity. Let us assume the standard 
collisionless Boltzmann equation: 

9t V 7 (34) 

-(V$- V„)/(f,tf,t) = 0, 

where, according to the Newtonian theory, only the po- 
tential $ is present. Considering the f(R) Poisson equa- 
tions, given by Eqs. (|32|) and |33|) . also the potential \P 
has to be considered so we obtain the coupled equations 



V 2 ($ + _>) - 2aV 4 ($ - = 16ttG / f(r, v, t)dv (35) 



V 2 ($-*) + 3aV 4 ($-*) = -8ttG / f(r, v, t)dv . (36) 

In the previous equations, we have replaced /"(0) with 
the greek letter a. It is important to stress that while in 
the standard theory, the distribution function f(f,v,t) is 
related only to the potential <!>, it is related to both $ and 
^ in the Newtonian limit coming from f(R) gravity. As 
in standard case, we consider small perturbation to the 
equilibrium and linearize the equations. After we write 
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equations in Fourier space so they became 

-iufi + v ■ (ikh) - (ik^) ' H - 0, (37) 

-fc 2 ($i + *i) - 2afc 4 ($i - *i) = 16ttG J fidv,(38) 

fc 2 ($i - *i) - 3afc 4 ($i - *i) = 8nG J f x dv. (39) 

Combining Eqs. (|5"5|) and (|3U)) . we obtain a relation be- 
tween $! and 

3 - Aak 2 ^ 
1 l-4afc 2 1 



inserting this relation in Eq. (|38j) and combining it with 
Eq. (|37|) . we obtain the dispersion relation 



1-4ttG 



1 - 4afc 2 
3afc 4 - A; 2 



dv 



dtf = 0. 



(40) 



If we assume, as in standard case, that fo is given by (flCl 
and fc = (fc, 0, 0), one can write 



1 + 



2V2ttG Po 1-iak 2 



3ak 4 - k 2 



kv x e 



kv x — ui 



-d,V x 



0. (41) 



By eliminating the higher-order terms (imposing a = 0), 
we obtain again the standard dispersion Eq. ^ . In order 
to compute the integral in the dispersion relation (|4T|). 
we consider the same approach used in the classical case, 
and finally we obtain: 



0.4 



0.2 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Figure 1. The bold line indicates the plot of the dispersion re- 
lation l|42[) in which we imposed the value for a given by (|43[) . 
The thin line indicates the plot of the standard dispersion 
equation [37l |. 



As shown in Fig. [TJ the effects of a different theory 
of gravity changes the limit of instability. The limit is 
higher than the classical case and the curve has a greater 
slope. This fact is important because the mass limit value 
of interstellar clouds decreases changing the initial con- 
ditions to start the collapse. 



l + g 



1 - 4afc 2 
3afc 4 - k 2 



1 - ^xe x2 (1 - erf [x]) = 0, (42) 



where x 



and Q = tt—. In order to evaluate 



V2ka <J 
Eq. (|42p comparing it with the classical one, given by 
Eq. §§§ , it is very useful to normalize the equation to the 
classical Jeans length showed in Eq. (|14p . by fixing the 
parameter of /(i?)-gravity, that is 



(43) 



This parameterization is correct because the dimension 
a (an inverse of squared length) allows us to parametrize 
as in standard case. Finally we write 



3fc 4 



fc 2 
k 2 



4P 



1 — y/nxe x (1 — erf[x]) 



(44) 

The function is plotted in Fig[U where Eq. (|42j) and the 
standard dispersion (37| are confronted in order to see 
the difference between f(R) and Newtonian gravity. 



V. THE JEANS MASS LIMIT IN /(i?)-GRAVITY 

A numerical estimation of the /(i?)-instability length 
in terms of the standard Newtonian one can be achieved. 
By solving numerically Eq. (|44|) with the condition uj = 
0, we obtain that the collapse occurs for 



1.2637kj. 



(45) 



However we can estimate also analytically the limit for 
the instability. In order to evaluate the Jeans mass limit 
in /(i?)-gravity, we set uj = in Eq. (|4"T|) and then 

(46) 



3cr 2 afc 4 - (WnGpoa + a 2 ) k 2 + inGpa = 0. 



It is worth stressing that the additional condition a < 
discriminates the class of viable f(R) models: in such 
a case we obtain stable cosmological solution and posi- 
tively defined massive states [3(J. In other words, this 
condition selects the physically viable models allowing to 
solve Eq. (|4"6"|) for real values of k. In particular, the above 
numerical solution can be recast as 



21 7T 



Gp 



(47) 
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The relation to the Newtonian value of the Jeans insta- 
bility is 



k 2 = 1 (3 + V2l) kj. (48) 
Now. we can define the new Jeans mass as: 



Mj = 6, 



(3 + \/2iy 



-M, 



(49) 



that is proportional to the standard Newtonian value. 
We will confront this specific solutions with some ob- 
served structures. 



V.l The Mj - T relation 

Star formation is one of the best settled problems 
of modern astrophysics. However, some shortcomings 
emerge as soon as one faces dynamics of diffuse gas evolv- 
ing into stars and star formation in galactic environment. 
One can deal with the star formation problem in two 
ways: i) we can take into account the formation of in- 
dividual stars and ii) we can discuss the formation of 
the whole star system starting from interstellar clouds 
To answer these problems it is very important 
to study the interstellar medium (ISM) and its proper- 
ties. The ISM physical conditions in the galaxies change 
in a very wide range, from hot X-ray emitting plasma to 
cold molecular gas, so it is very complicated to classify 
the ISM by its properties. However, we can distinguish, 
in the first approximation, between (40M43| : 

• Diffuse Hydrogen Clouds. The most powerful 
tool to measure the properties of these clouds is the 
21cm line emission of HI. They are cold clouds so 
the temperature is in the range 10 — 50 K, and their 
extension is up to 50-rT00 kpc from galactic center. 

• Diffuse Molecular Clouds are generally self- 
gravitating, magnetized, turbulent fluids systems, 
observed in sub-mm. The most of the molecular 
gas is H2, and the rest is CO. Here, the conditions 
are very similar to the HI clouds but in this case, 
the cloud can be more massive. They have, typi- 
cally, masses in the range 3 -7- 100 M@, temperature 
in 15 -T- 50K and particle density in (5 -=- 50) x 10 8 
m~ 3 . 

• Giant Molecular Clouds are very large com- 
plexes of particles (dust and gas), in which the 
range of the masses is typically 10 5 -7- 10 6 M Q but 
they are very cold. The temperature is ~ 15K. and 
the number of particles is (l-f-3) x 10 8 m -3 
However, there exist also small molecular clouds 
with masses M < 10 4 M Q ^. They are the best 
sites for star formation, despite the mechanism of 
formation does not recover the star formation rate 
that would be 250M Q yr~ 1 0]. 



• HII regions. They are ISM regions with temper- 
atures in the range 10 3 -7- 10 K, emitting primarily 
in the radio and IR regions. At low frequencies, ob- 
servations are associated to free-free electron tran- 
sition (thermal Bremsstrahlung) . Their densities 
range from over a million particles per cm 3 in the 
ultra-compact H II regions to only a few particles 
per cm 3 in the largest and most extended regions. 
This implies total masses between 10 2 and 10 5 M Q 
0. 

• Bok Globules are dark clouds of dense cosmic 
dust and gas in which star formation sometimes 
takes place. Bok globules are found within H II re- 
gions, and typically have a mass of about 2 to 50 
Mq contained within a region of about a light year. 



Using very general conditions [40 48], we want to show 
the difference in the Jeans mass value between standard 
and /(-R)-gravity. Let us take into account Eq. (fTB"]) and 
Eq.iPl): 



(50) 



Mr = -\ — 




in which po is the ISM density and a is the velocity dis- 
persion of particles due to the temperature. These two 
quantities are defined as 



Po = m H n H fj,, 



k B T 
m H 



where nu is the number 
fi is the mean molecular 
constant and ijih is the 
relations, we are able to 
interstellar clouds and to 
temperature. Results are 
have plotted the relation 
Any astrophysical system 
to a particular (M, T)-re, 
two theories for any self-g: 



of particles measured in m~ 3 , 
weight, ks is the Boltzmann 
proton mass. By using these 
compute the Jeans mass for 
plot its behavior against the 
shown in Tab|T] and FigH We 
for GR and for f(R) gravity, 
reported in Tab|I]is associated 
jion. Differences between the 
ravitating system are clear. 



Subject 


T n [l Mj Mj 

(K) (10 8 m- 3 ) (M ) (M e ) 


Diffuse Hydrogen Clouds 
Diffuse Molecular Clouds 
Giant Molecular Clouds 
Bok Globules 


50 5.0 1 795.13 559.68 
30 50 2 82.63 58.16 
15 1.0 2 206.58 145.41 
10 100 2 11.24 7.91 



Table I. Jeans masses derived from Eq. (|15|) (Newtonian 
gravity) and {49]) (/(ii)-gravity). 



By using Eq. (|49p and by referring to the catalog of molec- 
ular clouds in Roman-Duval J. et al. j49j, we have calcu- 
lated the Jeans mass in the Newtonian and f(R) cases. 
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20 30 40 50 60 
V[K] 



Figure 2. The Mj-T relation. Dashed-line indicates the New- 
tonian Jeans mass behavior with respect to the temperature. 
Continue-line indicates the same for /(i?)-gravity Jeans mass. 



Tab|TT] shows the results. In all cases we note a sub- 
stantial difference between the classical and f{R) value. 
In f(R) scenario, molecular clouds become sites where 
star formation is strongly supported and more efficient 
because in each of them the limit for the gravitational 
collapse is lower than the one in GR. 



VI. DISCUSSION AND CONCLUSIONS 

/(ii)-gravity is an approach aimed to address some 
shortcomings of modern cosmology just assuming exten- 
sions of GR without invoking the presence of dark in- 
gredients. In other words, dark energy and dark matter 
could be effects related to curvature further degrees of 
freedom instead of new fundamental particles. 

Here we have analyzed the Jeans instability mecha- 
nism, adopted for star formation, considering the New- 
tonian approximation of /(iJ)-gravity. The related 
Boltzmann-Vlasov system leads to modified Poisson 
equations depending on the /(i?)-model. In particular, 
considering Eqs.(|32|) and (|33|). it is possible to get a new 
dispersion relation (1421) where instability criterion results 
modified (see also [32|). The leading parameter is a, i.e. 
the second derivative of the specific /(i?)-model. Stan- 
dard Newtonian Jeans instability is immediately recov- 
ered for a = corresponding to the Hilbert-Einstein La- 
grangian of GR. In Fig. [TJ dispersion relations for New- 
tonian and a specific /(i?)-model are numerically com- 
pared. The modified characteristic length van be given 
in terms of the classical one. 

Both in the classical and in f(R) analysis, the system 
damps the perturbation. This damping is not associated 
to the collisions because we neglect them in our treat- 
ment, but it is linked to the so called Landau damping 
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Table II. We report the name, the particle number density 
and the excitation temperature of observed molecular clouds. 
For each system, we have calculated the value of Jeans mass 
in both Newtonian (GR) and f(R) gravity. The differences 
between the two approaches are significant pointing out that 
the star formation efficiency strictly depends on the adopted 
theory. This table is only a part of the catalog of molecular 
clouds reported in [49l |. 



A new condition for the gravitational instability is de- 
rived, showing unstable modes with faster growth rates. 
Finally we can observe the instability decrease in f(R)- 
gravity: such decrease is related to a larger Jeans length 
and then to a lower Jeans mass. We have also compared 
the behavior with the temperature of the Jeans mass for 
various types of interstellar molecular clouds (Fig. [2]) . In 
Tables[T]and|ll]we show the results given by this new limit 
of the Jeans mass for a sample of giant molecular clouds. 
In our model the limit (in unit of mass) to start the col- 
lapse of an interstellar cloud is lower than the classical 
one advantaging the structure formation. Real solutions 
for the Jean mass can be achieved only for a < and this 
result is in agreement with cosmology [30]. In particular, 
the condition a < is essentials to have a well-formulated 
and well-posed Cauchy problem in /(i?)-gravity [3(3]. Fi- 
nally, it is worth noticing that the Newtonian value is an 



s 



upper limit for the Jean mass coinciding with f(R) = R. 

This work is intended to indicate the possibility to deal 
with ISM collapsing clouds under different assumptions 
about gravity. It is important to stress that we fully re- 
cover the standard collapse mechanisms but we could also 
describe proto-stellar systems that escape the standard 
collapse model. On the other hand, this is the first step 
to study star formation in alternative theories of gravity 
(see also (33 - l36| ). From an observational point of view, 
reliable constraints can be achieved from a careful anal- 
ysis of the proto-stellar phase taking into account mag- 
netic fields, turbulence and collisions. Finally, addressing 



stellar systems by this approach could be an extremely 
important to test observationally /(i?)-gravity. 

Moreover, the approach developed in this work admits 
direct generalizations for other modified gravities, like 
non-local gravity, modified Gauss-Bonnet theory, string- 
inspired gravity, etc. In these cases, the constrained Pois- 
son equation may be even more complicated due to the 
presence of extra scalar(s) in non-local or string-inspired 
gravity. Developing further this approach gives, in gen- 
eral, the possibility to confront the observable dynamics 
of astrophysical objects (like stars) with predictions of 
alternative gravities. 



[1] T. Padmanabhan, Phys. Kept. 380, 235, (2003). 
[2] P.J.E. Peebles, B. Ratra, Rev. Mod. Phys. 75, 559, 
(2003). 

[3] V. Sahni., Led. Notes Phys. 653, 141, (2004). 
[4] S. Capozziello, M. De Laurentis, Phyis. Rept., 509, 167 
(2011). 

[5] S. Capozziello, M. Francaviglia, Gen. Rel. Grav. 40, 357, 
(2008). 

[6] S. Capozziello, M. De Laurentis, M. Francaviglia, S. Mer- 

cadante , Found, of Phys. 39, 1161, (2009). 
[7] S. Nojiri, SD. Odintsov, Int. J. Meth. Mod. Phys. 4, 115, 

(2007) . 

[8] S. Capozziello, M. De Laurentis, V. Faraoni, The Open 
Astr. Jour , 21874, (2009); Nojiri S., Odintsov S.D., 
Phys. Rept., 505, 59-144, (2011); S. Capozziello, V. 
Faraoni, Beyond Einstein Gravity: A Survey Of Grav- 
itational Theories For Cosmology And Astrophysics, 
Springer, New York, (2010); T. P. Sotiriou, V. Faraoni, 
Rev. Mod. Phys. 82, 451 (2010); S. Capozziello, M. De 
Laurentis, Invariance Principles and Extended Gravity: 
Theory and Probes, Nova Science Publishers, New York, 
(2011); T. Clifto n, P. G. Ferreira, A. Padilla, C. Skordis, 
larXiv:1106. 24761 [astro-ph.CO] (2011). 
[9] H.J. Schmidt, Lectures in Mathematical Cosmology, 
0407095, (2004). 

[10] S. Capozziello S., Int. J. Mod. Phys. D 11, 483, (2002). 

[11] S.M. Carroll, V. Duvvuri, M. Trodden, M.S. Turner, 
Phys. Rev. D, 70, 043528, (2004). 

[12] G. Cognola, et al., Phys. Rev.D 77, 046009, (2008). 

[13] G. Cognola, S. Zerbini, Int. J.Theor. Phys. 47, 3186, 

(2008) . 

[14] S. Nojiri, SD. Odintsov, Phys. Lett. B 576, 5, (2003); 

Phys. Rev. D 68, 12352, (2003). 
[15] R. Saffari, S. Rahvar, Phys. Rev. D, 77, 104028, (2008). 
[16] S. Nojiri, S. D. Odintsov, Phys. Rev. D, 77, 026007, 

(2008). 

[17] S. Capozziello, VF. Cardone, A. Troisi, Mon. Not. R. 

Astr., 375, 1423, (2007) 
[18] S. Capozziello, E. De Filippis, V. Salzano, Mon. Not. R. 

Astr. 394, 947, (2009). 
[19] Ch.G. Boehmer , T. Harko, F.S.N. Lobo , Astropart. 

Phys. 29, 386, (2008). 
[20] S. Mendoza, Y.M. Rosas-Guevara, Astron. Astroph, 472, 

367, (2007). 

[21] Y. Sobouti, Astron. Astroph. 464, 921, (2007). 

[22] S. Capozziello, S. Tsujikawa, Phys. Rev. D, 77, 10750, 



(2008). 

[23] S. Capozziello, M. De Laurentis, S. Nojiri, SD. Odintsov, 
Gen. Rel. Grav. 41, 2313, (2009); Physical Review D 79, 
124007, (2009); M. De Laurentis, S. Capozziello, Nuclear 
Physics B (Proceedings Supplements), 212, (2009); M. 
De Laurentis, S. Capozziello, S. Nojiri, SD. Odintsov, 
Annalen der Physik 19, 347, (2009); S. Capozziello, N. 
Carlevaro, M. De Laurentis, M. Lattanzi, G. Montani, 
larXiv:1104.2T69V l [astro-ph.CO] (2011); O.M. Lecian, 
G. Montani, Class. Quant. Grav. 26, 045014, (2009). 

[24] L. Pogosian, A. Silvestri, Phys. Rev. D 77 , 023503, 
(2008). 

[25] C. Bogdanos, S. Capozziello, M. De Laurentis, S. 

Nesseris, Astropart. Phys. 34, 4, (2009). 
[26] M. De Laurentis, S. Capozziello, Astropart. Phys. 35 257 

(2011). S. Capozziello, R. Cianci, M. De Laurentis, S. 

Vignolo S., Astroph. and Space Science 70, 341, (2010); 

S. Capozziello, M. De Laurentis, M. Francaviglia Astrop. 

Phys. 29,125, (2008); S. Capozziello, C. Corda, M. De 

Laurentis, Phys. Lett. B 669, 255, (2008). 
[27] A. Nishizawa et al., Phys. Rev. D 79, 082002, (2009). 
[28] K.S. Stelle, Gen. Rev. Grav. 9, 343, (1978). 
[29] CM. Will, Living Rev. in Relativity 9, 3, (2006). 
[30] S. Capozziello, S. Vignolo, Class. Quant. Grav. 26, 

175013, (2009). 

[31] JD. Barrow, A.C. Ottewill, J. Phys. A: Math. Gen., 16, 
2757, (1983). 

[32] S. Capozziello, M. De Laurentis, SD. Odintsov, A. Sta- 
bile, Phys. Rev. D, 83, 064004 (2011). 

[33] Ph. Chang, L. Hui, 2011, ApJ, 732, 25 

[34] A. Cooney, S. DeDeo, D. Psaltis, Phys. Rev. D, 82, 
064033 (2010). 

[35] A. Upadhye, W. Hu, Phys. Rev. D, 82, 064002 (2010). 
[36] C. A. Agon, J. F. Pedraza, J. Ramos-Caro, Phys. Rev.D 

83, 123007 (2011). 
[37] J. Binney, S. Tremaine, Galactic Dynamics, Princeton 

University Press, Princeton (1994). 
[38] L. D. Landau, E. M. Lifsits , Theoretical physics vol. 

II, Butterworth- Heinemann, Linacre House, Jordan Hill, 

Oxford (1987). 

[39] S.Weinberg "Gravitation and Cosmology", John Wiley & 
Sons, Inc., New York, (1972). 

[40] B.W. Carroll, D.A Ostlie., An Introduction to Modern 
Astrophysics, Pearson (2007). 

[41] R. Kippenhahn, A. Weigert, Stellar structures and evo- 
lution, Springer-Verlag, Berlin, (1990). 



9 



[42] M.A. Dopita, R.S. Sutherland , Astrophysics of the Dif- 
fuse Universe Springer, New York (2003). 

[43] H. Scheffler, H. Elsasser, Physics of the Galaxy and In- 
terstellar Matter Springer, New York (1982). 

[44] C. F. McKee The Origin of Stars and Planetary Systems, 
ed. C. J. Lada & N.D. Kylafis (Dordrecht: Kluwer), 29, 
(1999). 

[45] L. Blitz, Giant molecular clouds. In Protostars and Plan- 
ets III, eds. E.H. Levy and J.I. Lunine, (Tucson: Univ. 
of Arizona Press), pp. 125-161 (1993). 



[46] L. Blitz, L. Magnani, L. Mundy, Astrophys. J. Lett. , 
282, L9, (1984). 

[47] C. F. McKee, E.C. Ostriker, Annu. Rev. Astron. Astro- 
phys. 45, 565, (2007). 

[48] L.D. Anderson, T.M. Bania, J.M. Jackson et al. , ApJS 
181, 255A, (2009). 

[49] J. Roman-Duval, J.M. Jackson, M. Heyer, HJ. Rath- 
borne, R. Simon ApJ 723, 492, (2010). 



